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Abstract 


A summary  is  given  of  cur  present  state  of  knowledge  concerning 
optical  and  thermal  transitions  of  electrons  or  holes  trapped  in  solid 
Among  the  topics  covered  are  integrated  intensities  for  absorption, 
the  breadth  of  absorption  lines  and  the  rates  for  thermal  ionization 
and  capture.  The  relation  between  the  latter  and  photoconduct ive 
lifetimes  is  discussed  briefly. 
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I . Introduction 


The  scope  of  the  present  pax>er  will  be  limited  to  those  electronic 
transitions  in  solids  that  involve  lattice  vibrations  in  an  important 
way.^  The  simplest  category  is  that  of  radiation-Hess  transitions— 
the  thermal  ionization  of  a trapped  electron  (or  the  Inverse  process  of 
thermal  capture).  In  this  process  the  trapped  electron  through  its 
interaction  with  the  lattice  acquires  energy  by  absorbing  one  or  more 
phonons  and  is  released  into  the  conduction  band.  For  impurities  in 
germanium  or  silicon,  the  ionization  energy  of  .01  or  .0^  e.v.  is  less  . 
than  the  Debye  energy  k©  for  these  materials,  so  that,  in  principle,  a 
one-phonon  transition  can  take  place.  For  an  F center,  however,  the 
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ionization  energy  '•'ofAj  2 e.v.  In  KBr  is  many  times  larger  than  the 
energies  of  the  phonons  /\j  .03  6.v.,  so  that  thermal  ionization  can  only 
occur  as  a multi-phonon  process.  In  sections  2 and  3 explanation 
will  be  given  of  the  way  in  which  multi-phonon  transitions  are  mads 
possible. 

Although  optical  absorption  by  trapped  carriers  can  occur  without 
lattice  vibrations,  such  absorption  would  have  a line  spectrum.  The 
broadening  of  such  a line  into  a band  is  of  obvious  importance  for 
photoconductivity,  and  fur  this  lattice  vibrations  are  essential. 


An  attempt  will  be  mode  to  summarize  the  present  state  of  our 
knowledge  concerning  such  transitions.  References  will  be  given  to 
the  previous  literature  to  facilitate  the  reader  in  following  our 
discussion,  but  no  attempt  will  be  made  to  review  all  of  the  pr:i.io,:s 
literature,  and  significant  references  will  undoubtedly  be  omitted— 
with  apologies  to  the  authors  concerned. 

*'  R.  W.  Pohl,  Proc.  Pbys.  Soc.  42»  part)#  3 (1937)* 

3 B.  Szigeti,  Trans.  Far.  Soc.  155  (19**9)« 


u 

Mott  and  Gurney  using  the  classical  Pranck-Condon  principle 

suggested  that  the  broadening  should  increase  in  proportion  to  the 
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square  root  of  the  absolute  temperature.  Many  experiments'^  Indicate 

good  agreement  with  this  suggestion  at  high  temperatures.  At  low 

temperatures,  however,  the  line  width  did  not  approach  zero,  but  reduced 
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to  a constant  value.  For  KBr,  the  residual  line  width  of  about 

.2  e.v.  is  one  order  of  magnitude  larger  than  the  Debye  energy! 

12 

Schon  suggested  that  the  residual  line  width  was  produced  by  the 
zero-point  oscillations  of  the  lattice.  Calculations  based  on  a simple 
ccnfigurational  coordinate  model  (for  different  materials)  were  made  by 
Williams,'*'^  Klick,'*'^  and  Vlc.m'*’^  to  obtain  more  detailed  Information 
concerning  the  shape  of  the  bands  and  the  way  their  width  changed  with 
temperature . 


N.  F.  Mott  and  R.  W.  Gurney,  "Electronic  Processes  in.  Ionic  Crystals," 
Oxford  Press,  New  York,  19^0,  p.  117 ♦ 

^ C.  C.  Klick,  J.  Phys.  Chem.  776  (1953)* 

^ C.  C.  Klick,  J.  Opt.  Soc.  4l,  8l6  (1951). 

^ C.  C.  Klick  and  J.  H.  Schulnan,  J.  Opt.  Soc.  Am.  509  (1950). 

^ H.  Brinkman  and  C.  C.  Vlam,  Physica  650  (1949)* 

^ C.  C.  Vlam,  Physica  609  (1949). 

40  H.  P.  Hameka  and  C.C.  Vlam,  Physica  1£,  943  (1953). 

D.  D.  Johnson  and  F,  J.  Studer,  Phys.  Rev.  82,  970  (I951). 

M.  Schon,  Ana.  Physik  Series  6,  343  (1948). 

F.  E.  Williams,  J.  Chem.  Phys.  457  (l95l). 
l4 

F.  E.  Williams,  Phys.  Rev.  &2,  28l  (1951). 

F.  E.  Williams  and  H.  H.  Hibb,  Phys.  Rev.  84,  II81  (1951). 

C.  C.  Klick,  Phys.  Rev.  85,  154  (1952). 

C.C.  Vlam,  '‘The  Structure  of  the  Emission  Bands  of  Luminiscent  Solids 


Thesis,  University  of  Groningen,  July  1953. 


In  Section  3 we  shall  discuss  the  configurational  coordinate  model 

l6 

in  che  light  of  a generalized  Franck-Condon  principle.  In  particular, 

we  shall  discuss  the  hand  shape  and  width  at  low,  high  and  intermediate 

19 

temperatures.  We  shall  also  explain  the  observation  that  total 
intensities  of  absorption  (cr  emission)  are  essentially  independent 
of  temperature. 

ifi 

M.  Lax,  J.  Chen.  Phys.  1752  (1952).  Equations  quoted  from  this 
paper  will  be  prefixed  by  the  letter  I. 
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See  especially  reference  17,  figures  9 iL,  pp.  26-3?. 


2.  Optical  Absorption  by  Trapped  Electrons;  Integrated  Intensities 

An  electron  trapped  in  a crystal  will  have  a wave-function  whose 
extent  is  small  cooijared  to  the  wave-length  of  the  radiation  it  can 
absorb.  It  will  always  be  permissible  therefore  to  treat  radiation  or 
abstrption  by  making  a multipole  exj^ansion  of  the  electromagnetic  field. 
When  electric  dipole  radiation  does  not  vanish  because  cf  a selection 
rule,  it  will  make  the  dominant  contribution  and  all  higher  order  effects 
may  be  neglected. 

For  electrons  that  are  closely-bound  ("deep  traps"),  Including 

F centers,  the  electronic  energy  will  bo  large  c.nparpd  to  the  Debye 

energy  of  the  crystal.  The  electronic  frequencies  will  then  be  fast 

compared  to  the  highest  vibrational  frequencies  and  a Born-Oppenheimer 

treatment  of  the  electronic  notion  will  be  valid.  In  silicon  and 

germanium,  however,  the  (hydrogenlike)  impurities  have  ionization 

energies  cf  the  order  .05  e.v.  and  .01  e.v.  which  are  smaller  than  the 

Dehye  energy  and  the  Born-Oppenheimer  approximation  cannot  be  made 

without  further  analysis.  In  treating  the  broadening  of  impurity  levels 

in  semi-conductors,  it  was  found,  however,  that  the  important  vibrations 

are  those  whose  wave  lengths  are  ccmparahle  to  the  Bohr  radius,  more 

precisely,  whose  propagation  c»  nstant  = (2/a)  where  a is  the  Bohr 
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radius  in  the  crystal.  The  ratio  of  vibrational  energy  to  ground  state 
binding  energy  is  then  ho/Eg  ^ ^ (2/a)/Eg  where  is  the 

velocity  of  sound.  This  ratio  is  about  .1  or  .2  in  Ge  or  Si,  so  that  it 
is  still  a fairly  good  apprc-xlmation  tc.  employ  the  Born-Oppenheimer 
method . 


M.  Lax  and  E.  Bu~stein,  Broadening  of  Impurity  Levels  in  Silicon. 
To  he  submitted  to  the  Physical  Review. 


_4- 


An  analysis  of  the  absorption  and  enissicn  of  electric  dipole 

radiation  by  an  electron  trapped  in  a solid  has  already  been  given  using 

18 

the  Born-Oppenheiner  GX3prc,xir.iation.  The  principal  result  of  this 

id 

analysis  is  incorporated  in  the  starting  equation  l(2.l)  for  the 


absorption  cross-section: 

O-fv/)  - K \ F / ^ 


c 


_L  i^) 

ba 


(2.1) 


t 


The  factors  in  the  brackets  take  account  of  the  nodification  of  the 

electromagnetic  field  in  the  crystal:  n = optical  index  of  re±x-action, 

K = optical  dielectric  constant,  E^/E  = ratio  of  average  electric  field 

18 

at  the  center  tt.  its  average  value  in  the  crystal.  The  shape  of  the 

absorption  apectrun  is  contained  in  \)  ) where  the  subscripts  indicate 

that  a transition  has  taken  place  from  electronic  state  (b 

1 . V To. 

to  electronic  state  (with  any  accompanyiag  change  in 

the  vibrational  states).  Here  r is  the  electron  coordinate  and  x the 

set  of  nuclear  coordinates.  The  electric  moment  for  the  transition  is 


given  by 


(2.2) 
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The  integrated  absorption  from  1(4.3)  is  given  by: 

U'  l|»' 


IJ.y)d(hp)  ||^(x)ir^balx(cix 


(2.3) 


when  p^(x)  dy  is  the  (quantum-mechanical)  thermal  average 
probability  of  finding  the  nuclear  cov^rdinates  in  dx.  To  a zeroth 
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apprcximaticn,  M (x)  will  not  depend  significantly  on  x for  the 
ta 

noderate  displacements  that  occur  in  lattice  vibrations.  Then  (2.3) 


independent 


reduces  approximately  to  a constant  value 
cf  tenperature--in  agreement  with  the  widespread  cbservation  that  the 
integrated  absorption  is  essentially  constant  over  a wide  temperature 
range  even  though  the  breadth  and  shape  of  the  absorption  spectrum  may 
change  aijpreciably.  To  a first  approximation  C + Dx  and 


BO  that  the  effect  on  the  integrated  intensity  is  tc  produce  an 


increase  with  temperature  that  is  second  order  in  the  ratio  D/C.  (The 

primary  effect  of  such  a linear  term  is  to  alter  somewhat  the  shape 
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of  the  absorption  spectrum.  Dexter  has  found  that  the  usual  Gaucsian 
shape  can  be  altered  by  20  per  cent  or  more  on  the  wings.) 

The  integrated  absoirption  cross-sectinn  In  the  presence  of  vibration 
can  thc;refore  be  estimated  from  the  corresponding  cross-section  in  the 
non-vibrating  lattice.  We  can  avoid  statements  tliat  depend  on  an 
explicit  choice  of  electronic  wave-functions  by  introducing  the  oscillatoa 
strength 

L-  - E J C'/3)  \ It  (...) 


The  electronic  (rigid-lattice)  cross-section  to  all  possible  final  states 
b can  then  he  written 


D.  L.  Dexter,  On  the  Shapes  of  Absorption  and  Emissicn  Lines  of 
Impurities  in  Solids,  Office  of  Scientific  Research  Technical  Note 

54-192= 

<■ 

-o- 


where  1 {\j)  — ^ 

' \)70,  ^ 

jrU)ci(KD)-X  ibO- 

A,>(X 


Each  tenn  in  (2.6)  represents  a sharp  line  because  of  the  neglect 

of  all  sources  of  broadening,  however,  the  corresponding  torn  in  (2.7) 

yields  the  integrated  result  for  that  line,  which  to  a first  approxioatior 

is  valid  even  in  the  presence  of  broadening. 

Note  that  □ cancels  out  of  the  codbiued  equations  (2.U^)  and  (2.5)> 

and  can  therefore  be  chosen  sornewhat  arbitrarily.  For  a tiighiy  localised 

electron  e.g.  in  an  F center,  the  first  allowed  transition  (e.g.  Is  to 

« 

2p)  will  have  an  oscillator  strength  of  the  order  unity  (using  in  m) 
and  higher  transitions,  e.g.  Is  to  np  will  be  reduced  roughly  by  a 
factor  n-'  because  cf  the  reduction  in  overlap.  Thus  (2.5)  and  (2.?) 
with  these  assumed  oscillator  strengths  provide  a quick  rough  estimate 
of  the  integrated  cross-sections  for  the  discrete  transitions.  Direct 
calculations  of  the  oscillator  strengths  are  necessary,  for  a more 
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accurate  analysis.  Such  calculations  have  been  carried  out  by  Dexter 
for  the  F-center  case. 

For  electrons  whose  wave  x'unctioas  are  sufficiently  slowly  varying, 
the  Wannier-Slater-Janies  treatment  of  electrons  in  perturbed  periodic 
lattices  becomes  valid.  It  then  becomes  possible  to  represent  the  wave- 
function  cf  a trapped  electron  as  a linear  combination  of  orthogonalized 
atomic  (Wannier)  orbitals  belonging  to  one  band  (if  no  degeneracy  exists). 
The  coefficients  in  this  expansion  are  a slowly  varying  function  of  their 
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D.  L.  Dexter,  Phys.  Rev.  435  (1951)- 
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argument— the  cell  position,  and  when  replaced  by  a single  continuous 
function  describe  the  smooth  envelope  of  the  complete  wave  function 
whose  rapid  variations  are  contained  within  the  Wannier  orbitals. 

This  smooth  envelope  function  may  be  regarded  as  a macroscopic  Schrodingei 
wave  function  since  it  bears  the  sane  relation  to  the  complete  wave-, 
function  as  a macroscopic  electromagnetic  field  does  to  the  microscopic 
field.  Indeed,  the  macroscopic  Schrodlnger  function  obeys  a hydrogen- 
like wave  function  with  an  effective  mass  different  from  the  vacuum 
mass  just  as  a macroscopic  electric  f iel'?  obeys  a Maxwell  equation 
containing  the  medium  dielectric  constant. 

A slight  extension  of  the  Wannier -Slater -Janes  treatment  can  be 
used  to  show  that  for  electrons  with  large  orbits,  the  electric  dipole 
matrix  element  can  be  calculated  to  a good  approximation  directly  from 
the  macroscopic  Schrodlnger  wave  functions  ignoring  the  periodic 
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substructure.  A consequence  cf  this  result  is  that  the  usual  sun 
rule  for  oscillator  strengths 


au  b 


(2.8) 


is  valid  providing  that  n in  (2.U)  is  the  effective  mass  in  the  crystal. 

If  ellipsoidal  energy  surfaces  are  present  then  (2.8)  remains  valid 

* 2t 

provided  a Is  interpreted  as  the  harmonic  mean  of  the  three  masses: 


l/m*  = l/mi  + l/nig  + 


(2.9) 
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See  for  example  how  the  matrix  elements  of  the  deformation  potential 
are  calculated  in  ref.  I8. 

2k 

A proof  of  this  statement  will  be  included  in  a paper  by  E.  Burstein 
and  the  author  cn  optical  absorption  by  impurity  levels  in  silicon. 


I 


If  the  energy  surfaces  are  spherical  and  the  potential  hydrogen- like, 
then  the  oscillator  strengths  are  the  sane  as  in  the  original  hydrogen 
aton — the  inl'lusnce  of  a change  in  effective  mass  and  dielectric  constant 
cancels  cut  of  the  oscillator  strength  fomula  (2.i|)  since  E.  -E„  is 

O 


* . 2 


proportional  to  n /K  and  r^^  is  proportional  to  the  modified  Bohr  radius, 

/ “1 

or  to  {Kin  ) 


The  absorption  in  the  continuous  spectrum  also  bears  a simple 
relation  to  that  in  a hydrogen  atom.  The  contribution  of  the  continuous 
'"(u)  can  be  written: 


spectrum  to 


I 


_./  I 


\ 2^  I ] ' ji//p'^Tn\  CulCXOi  nvj  (p  10) 


t-ONT' 


V'!  3 i U u i r-  2TT  m'J 


»herc  V / t “ 1 J ^ 


and  Eg  is 


the  ionization  energy  of  the  ground  (is)  state.  The  integral 

pj=  .U30  represents  the  fraction 
of  the  absorption  iiT  the  continuum  as  compared  to  the  total  absorption. 

The  contribution  to  this  integral  from  any  region,  ’Measured  in  units  of 
Eg,  is  identical  to  the  corresponding  contribution  for  the  hydrogen  atom. 
However,  ( XJ  ) itself  according  to  (2.5)  and  (2.10)  is 
proportional  to  (n/K)(n  Eg)  or  IT'  /(r.  ) . 

For  purposes  of  comparison  with  experiment,  we  summarize  the  most 
pertinent  hydrogen-like  oscillator  strengths  from  Is  to  np  and  to  the 
continuum:*^^ 


%,1S  = f3p,ls  - -079,  fUp.ia  = -029,  = -kSP  (2.U) 


H.  Bethe,  Handbuch  der  Phys.  2U/l,  Kap.  3>  Ziff.  4l,  Julius  Springer, 
Berlin  1933- 
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and  the  integrated  ahsorption  croBS-section:  ^ 

r , ^ 4r'  o = ^TT’'t^e 

j O'(u)  OVi'W  VE  / K,  frrf  mi  c 


(u)  d(nu)  - (^) 


(2.12) 


vith  — d. — " ^ — 1*098  X 10”^°  e.v.  (cm)  . 

figure  1 and  2 show  the  optical  absorption  by  Boron-doped  silicon. 
Three  absorption  peaxs  can  be  observed  corresponding  to  the  ls-2p,  ls-3p 
and  Is-Up  transitions  in  addition  to  the  continuum  absorption  above 
about  0.046  e.v.  The  most  striking  discrepancy  betv'een  theory  and 
experiment  is  that  the  discrete  absorption  peaks  get  about  10  per  cent 
of  the  total  absorption  rather  than  about  the  50  per  cent  to  be  expected 
from  the  hydrogen  model.  Secondly,  the  strongest  discrete  transition 
is  the  1s-3p  transition  rather  than  the  ls-2p  transition.  These 
discrepancies  are  related  to  the  fact  that  the  valence  band  of  silicon 
contains  three  degenerate  bands  (six  if  spin  is  included),  and  spin-orbit 
coupling  only  partially  removes  the  degeneracy.  Thus  the  wave  functions 
describing  an  impurity  trapped  hole  involves  a linear  combination  of 
Wannier  functions  from  several  bands,  and  the  elementary  treatment  given 
above  is  not  valid. 

The  integrated  absorption  should,  however,  be  less  sensitive  to 
details  of  the  wave-functions  and  better  agreement  should  be  obtained 
with  the  theoretical  expression  (2.12).  The  experimental  integrated 


absorption  coefficient 


Joi(V)  d(hl>) 


is  0.90  (e.v.)cm  . 


This  may  be  converted  to  an  integrated  cross-section  by  dividing  by  the 
concentration  of  Boron  centers.  To  estimate  the  latter,  we  note  that 
the  specimen  in  question  has  a room  temperature  resistivity  of  1 ohm-cm. 
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Assuming  a mobility  of  350  cn'^/volt-sec,  this  leads  to  a room 
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temperature  hole  concentration  of  1.8  x 10  holes/cm'^.  Since  all  Boron 

2 

centers  would  be  ionized  at  room  temperature,  we  nay  take  1.8  x 10^°/cm 
for  the  concentration  of  Boron  centers--assuming  the  concentration  of 
other  centers  to  be  negligible.  Thus  we  can  conclude  that 

CH  0.5  X 10"^^  e.v.(cm)^  (2.13) 

with  an  accuracy  no  better  than  30  per  cent. 


The  theoretical  integrated  absorption  cross-section  can  be  estimated 

2 

from  (2.12)  assuming  (l)  the  effective  field  correction 

negligible  for  a center  as  large  as  the  one  in  question  (2)  the 

effective  mass  may  be  determined  from  the  ionization  energy  Eg  Or  .046  e.v. 

^ . .-1/2^  . .-1/2 


to  be  m Cir:  0.45  m (3)  n/k (k)  ' (12)  ' . The  final  result  of 

-16  ^ 2 

0.7  X 10  ev.  (cm)  is  In  reasonable  agreement  with  experiment. 
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This  mobility  has  been  found  reasonably  representative  at  Bell 
Laboratories  for  estimating  carrier  concentrations  in  a group  of 
specimens  from  which  the  sample  in  question  was  taken. 


3.  Absorption  by  Traijped  Electrons:  Broadening 

There  are  two  distinct  stages  in  the  calculation  of  the  shape  of  an 

absorption  curve  broadened  by  lattice  vibrations.  TSe  first  is  a 

determination  of  the  electrtnic  energy  as  a function  of  the  i>ertinent 

nuclear  coordinates  using  either  one  configurational  coordinate  or 

many.  The  second  is  the  calculation  of  the  spectrum  from  the  known 

energy  v«.  configurational  coordinate  curves  using  classical,  semi- 

classical,  or  qu'^ntum  mechanical  methods. 

Simple  energy  level  diagrams  against  a single  fictitious  coc>rdinate 

have  long  been  used  to  give  a qualitative  understanding  of  absorption, 

27 

luminescence,  non-radiatlve  transitions,  etc.  The  first  serious 

attempt  to  calculate  accurately  the  energy  against  a (single)  real 

28 

coordinate  was  made  by  Williams  fc-i  thallium  activated  KCl.  Williams 
used  as  his  coordinate  the  radial  displacement  of  the  six  nearest 
neighbors.  His  configurational  diagram  is  shown  in  Figure  3«  The 
remarkable  agreement  between  his  calculated  spectrum  and  the  observed 
spect— am  is  shown  in  Figure  4.  Williams  results  will  be  analyzed  in 
more  detail  later. 

27 

See  for  example,  H.  W.  Leverenz,  Luminescence  of  Solids,  John  Wiley 

and  Sons,  New  York  1950,  p.  132,  or  Mott  and  Gurney,  ref.  4. 

28 

P.  £.  Williams,  J.  Phys.  Chem.  78O  (1953)*  This  paper  contains 
the  most  accurate  calculation  of  the  configurational  diagram.  More 
detailed  discussion  is  given  in  the  earlier  ixxpers,  refs.  I3  - 15* 
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Klick  and  Vlan  have  eseentially  reversed  the  Williams  procedure, 
hy  deriving  configurational  coordinate  curves  from  the  experimental  data. 
They  assume  that  an  adequate  description  can  he  given  in  terras  of  one 
coordinate,  and  that  the  energies  of  the  ground  and  excited  states  nay 
he  approximated  hy  quadratics  in  the  coordiiiate.  Nn  attempt,  however, 
is  nade  to  assign  a real  significance  to  the  "fictitious"  coordinate,  or 
to  compare  the  configuration  cum^^es  with  calculations  from  fundamental 
principles.  Nevertheless,  a significant  contrihuticn  has  bc-en  made  hy 
this  reduction  of  the  experimental  data  to  a more  basic  form.  The 
configuration  curves  permit  the  observed  absorption  (anu  emission)  to  be 
qualitatively  understood,  and  provide  an  intermediate  link  between  experi- 
ment and  a fundamental  theory.  For  the  many  cases  in  which  the  basic 
mechanism  for  broadening  is  not  rnderstood,  or  too  difficult  to  treat, 
the  Klick  and  Vlam  procedure  is  the  best  that  can  be  done. 

29 

Two  courageous  physicists  Huang  and  Rhys  in  treating  the  F center 
problem,  attempted  the  double  feat  of  (l)  calculating  the  F center  energy 
ns  a function  of  all  the  normal  cotrdinates  in  the  lattice  (2)  making 
an  all  quantum  mechanical  treatment  of  the  many  body  problem.  As  we 
shall  see,  they  were  mere  successful  in  the  second  or  raethodologicol 
aspect  of  their  work. 

Huang  and  Rhys  pointed  out  that  the  energy  change  of  the  F center 
electron  would  he  produced  primarily  by  the  long  range  electrostatic  forces 
present  in  a polar  lattice.  Furthermore  they  noted  that  large  long- 
range  fields  are  produced  primarily  by  the  long-wave  longitudinal 
optical  modes.  They  therefore  neglected  interactions  with  all  other 

2Q  ^ 

K.  Huang  and  A.  Rhys,  Proc.  Roy.  Soc.  (London)  A20k,  ko6-423  (1950). 
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inodes,  and  treated  the  long-vave  longitudinal  optical  inodes  by  a 

30 

dielectric  continuun  model  earlier  discussed  by  Huang.  This  continuum 

model  smooths  the  cell  polarization  into  a continuous  function — just 

as  the  deformation  potential  method  smooths  the  displacements  in  a 

long  wave  acoustic  mode.^^ 

29 

Huang  and  Rhys  describe  the  electron-lattice  interaction  as  the 
interaction  of  the  continuum  polarization  with  the  electric  field  produced 
by  the  static  charge  distribution  associated  with  the  F center  electron. 
Since  the  electric  field  decreases  very  slowly  with  distance,  Huang  and 
Rhys  expect  a long  range  induced  polarization,  which  when  translated 
into  normal  coordinates  emphasizes  the  long  wave  length  modes. 
Unfortunately,  however,  Huang  and  Rhys  neglected  the  fact  that  the  F 
center  is  neutral.  Stated  differently,  the  only  polarization  that  we 
need  consider  is  the  change  compared  to  the  perfect  lattice  case,  and 
therefore  the  charge  distribution  to  be  used  is  that  of  the  F center 
electron  minus  that  of  the  negative  ion  it  replaces.  Such  a neutral 
distribution,  however,  produces  only  short  range  effects  and  invalidates 
the  continuum  approximation  made  by  Huang  and  Rhys. 

One  might  anticipate  that  the  Huang  Rhys  procedure  of  using  the 
static  unneutralized  charge  distribution  of  the  electron  would  yield  too 
high  an  answer  for  the  broadening.  Actually  their  parameter  S which 
measures  the  broadening  is  calculated  to  be  3*6  whereas  a cemparison  with 


K.  Huang,  Report  L/T  239  of  British  Electrical  and  Allied  Industries 
Research  Association.  See  also  H.  B.  Callcn,  Phys.  Rev.  JC,  139^ 
(19^9)  and  II.  Frohlich,  Free.  Roy.  Soc.,  London  AI60,  23O  (1937)* 

J.  Bardeen  and  W.  Shockley,  Phys.  Rev.  72  (1950). 
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experinent  yields  S 22.  Two  points  are  involved  in  explaining 
their  small  result.  The  first  is  that  the  error  involved  in  using  a 
non-neutral  charge  distribution  cancels  out  of  their  final  calculations 
since  only  the  difference  between  the  electron- lattice  interactions  in 
the  ground  and  excited  states  is  of  importance.  The  second  is  that 
their  continuum  treatment  grcissly  underestimates  the  contributions  of 
short  range  interactions. 

To  explain  the  latter  statement,  we  note  that  the  electron-lattice 

interaction  used  by  Huang  and  Rhys  is  linear  in  the  normal  coordinates 

and  its  net  effect  is  to  shift  the  equilibrium  positions  of  all  the  ions. 

Hc.wever,  the  major  shift  will  take  place  for  the  nearest  neighbor  ions. 

and  such  a shift  cannot  be  calculated  from  a continuum  mc^el  but  must  be 

based  on  a detailed  analysis  of  all  the  local  forces  similar  tc  the  cue 

carried  out  by  Williams  fc-r  KCltTl. 

32 

Burstein  and  Oberly  have  suggested  a model  for  the  F center 

based  on  treating  the  trapped  electron  as  a particle  in  a box.  This 

33 

model  correlates  with  the  empirical  rule  suggested  by  Mollwo 

'Z.  _1  2 

y)  CX  '2:1.  sec  cm  (3-1) 

rrvKTjC 

relating  the  maximum  absorption  frequency  ) J to  the  inter-ion 

^max 

distance  a for  all  the  alkali  halides.  The  particle  in  a box  model 
emphasizes  the  incompressibility  of  the  alkali  ions  that  form  the  box 
walls — i.o.  it  emphasizes  short  range  forces  of  a sort  omitted  in  the 
Huang  Rhys  treatment.  A closely  related  molecular  model  has  been  proposed 

•32 

E.  Burstein  and  J.  J.  Oberly,  Optical  Properties  of  F Centers  at 
Liquid  Helium  Temperatures.  Proceedings  of  Conference  on  Low 
Temperature  Physics,  1951,  Natiofial  Bureau  of  Standards  Circular  519. 
33  E.  Mollwo,  Z.  Phys.  §3,  56  (1933). 
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l.r 


3U 

by  T.  Innui  and  Y.  Venura,  and  a detailed  molecular  calculation  was 
35 

started  by  Muto.  However,  an  adequate  calculation  of  the  configuration 
energy  of  the  F center  does  not  yet  exist. 

In  the  remainder  of  this  section,  we  shall  assume  that  the  energy 
of  the  trapped  electron  as  a function  of  the  various  configurational 
coordinates  is  known,  and  shall  discuss  the  methods  available  for 
computing  the  shape  of  the  absorption  spectrum.  The  Franck-Condon 
principle  states  that  we  any  regard  the  nuclei  as  standing  still  while 
the  electron  ciicingea  its  state  due  to  absorbing  the  photon.  Thus  the 
strength  of  absorption  can  be  computed  for  a given  configuration  and  the 
result  averaged  over  all  possible  configurations.  The  shape  factor 


X X)  in  (2.1) 


can  then  be  written 


(3.2) 


I \\XX 


(3.3) 


The  squared  matrix  element  1 Mb(o(^)l  represents  the  strength 
of  the  absorption.  F^(x)  dx  is  the  normalized  probability  of  finding  the 
nuclei  between  x and  x + dx  (here  x may  represent  a set  of  nuclear 
coordinates)  and  the  presence  cf  the  delta  function  factor  insures  that 
only  those  configurations  whose  energy  difference  d^l(x)  = V\  V will 
contribute  to  tl;e  spectrum  at  K u . A "derivation"  and  detailed 

16 

analysis  of  the  Franck-Condon  principle  has  been  given  by  the  author. 


T.  Innui  and  Y.  Uemura;  Prog.  Theor.  Phys.  395  (1950) 
T.  Muto,  Prog.  Thecr.  Phys.  U,  l8l  (19^9). 


We  shall  refer  to  the  Franck-Coadon  principle  as  classical  if 


classical  statistics  are  used  in  obtaining  the  probability  distribution: 


Pa(x)ciassical  “ 


Q Eo.U'i  At] 


h) 


and  as  seni-classical  if  the  quantun  mechanical  distribution  is  used: 


P M [■ 

■tc^yaATl'l  Xa.v^-(X)f 

' V’  r 

/m  1 

Eoym  / 

1 

(3.5) 


« . 4-t. 

where  is  the  m"*'  vibrational  wave  function  when  the  electron  Is 

in  state  a,  and  is  the  total  electron-vibrational  energy  in  this 
state.  For  sufficiently  high  temperatures,  the  quantum  mechanical 
distribution  (3*5)  can  be  shown  to  reduce  to  its  classical  value  (3.4). 
Such  a classical  approximation  will  be  valid  providing  kT  u) 

where  <a  is  the  frequency  of  the  normal  mode  or  modes  important  in  the 
broadening.  We  shall  see,  however,  that  the  semi-classical  raethcd  is 
so  easy  to  apply  that  little  is  gained  by  making  the  classical 
approximation . 

We  shall  illustrate  our  remarks  by  considering  the  case  where  x is 
a single  variable,  (For  the  appr'?priate  generalization  to  the  many  body 
problem,  the  reader  is  referred  to  reference  l8).  Equation  (3*5)  can 
then  be  reduced  to 


I (3.6) 
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and  if  the  states  X^„(x)  can  be  approxinated  by  states  of  a harnonic 


oscillator  of  frequency  co  and  nass  M then 
where  ~ 


and 


-1 


(3  7) 

(3.8) 

(3.9) 


The  conbination  of  (3.6)  and  (3*7)  yields  the  cciaplcte  spectral  distribu- 
tion. To  be  even  more  specific,  I shall  assume  that  the  excited  as  well 


as  ground  state  configurations  can  be  approxinated  by  parabolas: 

.a. 


AE(x)=A%-^6  A + ki), 


bO- 


(3.10) 


(3.11) 


where  ^ l^bo.  is  the  energy  difference  when  x = 0 i.e.  from  the 
equilibrium  point  of  the  lower  state  and  K - B is  the  stiffness  of  the 
excited  state.  (Usually  B ^ O). 

The  spectrum  (3-6)  aside  from  a factor  I (X)l  can  now  be 


written 

( 


O 


( v)  =[z 7T<xtf  ^ APj  Vft  1 ^ 

+ +r/' 


where  the  two  roots  for  x are  x=(.^/B)(l+p  ) and 
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If  the  change  in  stiffness  B is  sufficiently  small,  or  the  temperature 
sufficiently  low  tliat  then  it  is  permissible  to  drop 

the  second,  exponential  and  expand  with  the  result  that  the  spectrum 


becoines  Gaussian:  ^ ?*  / “1 

G(v)cr[rw<x^3^/\  M (3. lit) 


with  Wi'xS?'  ' ^ 


(3.15) 


and  a mean  energy  independent  of  temperature. 


A detailed  quariLuiu  L'leclianical  analysis  of  the  opeetx-ut;i  showed 
that  the  semi-classical  spectrum  (3.12)  has  the  correct  first  and  second 
moment ; 


<hy>  = kv  =<aE(x)>=  hv^-4:6<x^? 


<(Kp  - hvf  > = A<x^>+  X 


(Note  the  temperature  dependence  of  the  mean  energy  (3.1b)).  However, 
the  shape  of  the  spectrum  (3»12)  depends  on  the  validity  of  the  semi- 
classical  approximation,  which  we  shall  now  examine,  by  considering  the 
extent  to  which  the  third  moment  is  adequately  given  by  the  semi -classical 
formula 


The  difference  between  the  quantum  mechanical  and  semi -classical  third 

18 

moments  is  found  from  1(4,5)  and  (4.8)  to  be: 


(3.19) 
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where  T - P^/2M  is  the  kinetic  energy,  (T,i^!E)  is  a connutator,  and  the 
symbol  A'u-  means  take  an  average  of  the  enclosed  operator  in  the  state 
X^^(x)  and  then  take  a Boltznann  average  over  the  states  in.  It  is  not 
fair  to  coirpare  the  error  (3*19)  with  (3*18)  since  (3.18)  vanishes  if 
E = Oo  An  estimate  of  the  accuracy  involved  can  be  obtained  by  comparing 

(3.19)  with 

This  criterion  was  stated  in  1(4.9).  Indeed,  the  ratio  of  the  third  moment 
to  the  three  halves  power  of  the  second  moment  is  called  the  skewness  and 
denoted  ^ . (c.f.  Eq.  I 8.l8).  Our  criterion  for  validity  of  the 

semi-classical  approxiniation  is  M ?^S.C. 

If  this  condition  is  otc  yed  the-n  (3.12)  will  be  valid.  The  spectrum 
will  however  reduce  to  a Gaussian  even  in  the  semi -classical  case  only 
if  the  skewness  is  then  sufficiently  small:  25.C.«  1 

Evaluating  (3-19)  and  (3.18)  and  dividing  by  the  three-halves  power  of 
(3.17)  we  obtain  respectively: 

■ ^S.c. 


^sc.  - 


(3.21) 


2 — 2 

where  \ = (16/17)  (K/B)  /(l  + 2 /W  ) 


2,.  _ o . 2, 


(3.22) 

and  C<  = ^ a‘'/(17  -B‘^  <x‘‘>  ) (3.23) 

The  semi -classical  criterion  -1,,«  1 can  be  simplified 

by  considering  separately  the  cases  in  which  \ is  greater  or  less  than 
unity.  If  \>1,  then  c<>>  \ or  A^(2 /T\  +l)/(2  1\  Ito^)  » 1 (3.24) 

In  this  case,  it  is  clear  that  we  must  also  have  oC.  1,  so  that 

using  (3-21)  the  si>ectrum  is  also  Gaussian.  Comiparing  with  (3.8)  and 
(3.17);  the  condition  (3-24)  can  bp  written  in  the  more  perspicuoqs:' form. 
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(3.25) 


If  \>  1,  and  C^'>  1 it  fcllcwa,  ccnparing  (3.25)  and  (3.17)  that  the 
ohserved  broadenine  will  be  large  conpared  to  the  phcnon  energy  of  the 
ground  state,  (in  the  absence  of  other  information  the  Debye  energy 
k0  can  be  used  as  an  estimate  for  *1^  to) . This  latter  condition  is 
necessary  but  not  sufficient  for  the  validity  of  the  semi-classical 
method.  Even  if  \ 1,  if  ^ 1,  the  condition  (3«25)  is  sufficient 
to  insure  the  serai -classical  and  gaussian  nature  of  the  spectrum. 
However,  if  1,  the  semi-classical  condition  (3.20)  merely  requires 

- "^‘^3 

OC,  > ^ when  \ < 1 (3.26) 

1/3 

so  that  if  1 o(^  3>  \ a ser.il -classical  nen-gauesian  spectrum 

can  be  observed.  Under  these  circumstances  the  condition  (3.26)  can  be 
rewritten  in  the  form:  _ . / 

KUVo( 


(3.27) 

(3.28) 


» (txWo()[i  + 

which  is  slightly  more  stringent  than  (3-25)  since  1.  Therefore 

if  (3.25)  is  disobeyed,  a fully  quantum  mechanical  treatment  is  necessary. 

We  shall  illustrate  these  criteria  by  applying  them  to  the  problem 

of  KC1:T1  discussed  by  Williams.  From  his  configuration  diagram 

(Figure  2)  wc  may  deduce  that  A = 2.393  e.v./A,  A = 8.OI6  e.v./A'^, 

K = 16.85  e.v./A^,  ^ cu  = .0166  e.v.  using  his  mass  of  M = 6 
o I 22 

+ (.1264)  = 4.239  X 10  grans.  At  aosolute  zero,  we 

<x"f  2.2  2x|5^A 

ACO-  I U K’Xn  b')  - 1 U 

0^(0)  = 


have 


(3.29) 
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Thus  we  have  \ 'y  1 and  o<>»  \ > 1 so  that  the  spectrum  will  be 

both  semi-classical  and  Gaussian.  In  fact  k-l 
ana  ^^-2=  ( .02  so  that  the  skewness  would  be 

quite  small.  These  results  could  have  been  deduced  directly  from  the 
experimental  data  (assuming  X 1)  by  noting  that  the  r.m.s.  width 

i, 

<(bv  - .125  e.v./(2.346)  Ci=  .0533 

(where  .125  e.v.  is  the  full  width  at  half -power).  This  width  is 
larger  than  the  Debye  energy,  and  certainly  larger  than  o)  ^ .0166  e.v 


in  fact  ^ 


( .OI06/.O533)  — .1  in  agreement 


with  the  previously  calculated  value. 

2 

At  higher  temperatures  X.(T)  =»\(  O )/(l  + 2/v\)  , 

0<(T)  = o((  O ) / (1  + 2 7^)  so  that  \(T)/  decreases 

with  increasing  temperature,  and  the  semi -classical  criterion  remains 
valid  in  William's  case.  The  spectrum  will  remain  Gaussian  as  long  as 
O^T)  = (42.6)/(l  +•  2 /y\  ) > '^  1,  i.e.  for  all  measurable  temperatures 
As  a contrast  to  the  William's  case,  we  shall  con«»ider  the  configura- 
tional diagram  deduced  by  C.  C.  Klick  and  J.  A.  Schulman  for  Zn2SiO|^:Mn. 

-4  4 2 

They  find  A — -.9  x 10  ergs/cra,  B = 31  x 10  ergs/cra  , 

K = 46  X 10^  ergs/cm^,  M = 4m  ^ 1.057  x 10”^^g.,  “K  o>  ^6.952  x 
-14 

10  ergs  ^ .0434  e.v.  The  observed  broadening  at  low  temperatures 

is  of  the  order  of  .1  e.v.  so  that  one  might  hope  for  the  validity  of 
a semi -classical  treatment.  However,  the  zero  pt  displacement 
( "^(jo/2k)^^^0=:  2.67  X 10’^^  and  A(  "K  to/2k)^^^^  2.4  x 10  ergs 

which  is  only  .35  ( "K  <u)  so  that  condition  (3.25)  is  violated.  One 
finds  further  that  \(  O ) ~ 2,  and  0((  O ) ~ .27.  With  these 
parameters,  the  semiclassical  requirement  cannot  be  satisfied  at  any 
temperature.'  The  basic  reason,  of  course,  is  that  the  two  parabolas 
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of  the  configuration  diagram  have  appreciably  different  stiffnesses, 
yet  their  minima  are  quite  close  together  so  that  (a)  the  broadening 
is  produced  priaailly  from  the  difference  in  curvature  (b)  matrix  elements 
are  only  appreciably  to  low-lying  vibrational  states  and  for  these 
quantum  effects  are  appreciable. 

A comment  may  be  added  here  about  the  P center.  The  Huang-Rhys 
parameter  S is  simply  A /(  find  since  S is  given  the 

empirical  value  of  22,  condition  (3.25)  is  easily  satisfied.  Since  we 
may  anticipate  that  K/B  and  hence  \(o)  will  be  of  the  order  upity, 

(3.25)  is  a sufficient  condition  for  the  validity  the  semi  — C xclSolCaX 
method.  We  therefore  expect  that  the  F center  absorption  can  be 
adequately  represented  by  a Gaussian  at  low  temperatures,  and  by  (3-12) 
at  high  temperatures,  except  for  the  influence  of  cubic  terms  in  the 
latter  case. 

33 

Early  measurements  of  P center  absorption  due  to  Mollwo  are  shown 
in  Pigure  5*  More  careful  measurements  have  been  made  recently  by 

36 

Hesketh  and  Schneider  in  an  attempt  to  analyze  the  shape  of  the  F band. 
Figure  6 shows  a comparison  of  their  experimental  results  with  a 

37 

Lorentzian  shape— such  as  had  previously  been  assumed  by  Snakula  . and 
others  to  be  valid  on  the  basis  of  dispersion  theory.  The  agreement  is 
poor- -since  the  fundamental  mechanism  for  broadening  is  not  the  usual 
one  of  dispersion  theory.  Figure  7 shews  a Gaussian  plot  by  Hesketh 
and  Schneider  '.f  their  data.  A straight  line  of  slope  /y\.  indicates 

36 

R.  V.  Hesketh  and  E.  E.  Schneider,  Phys.  Rev.  95,  8372  (195^). 

-?7 

A.  Smakula,  Z.  Physik  603  (1930). 


-23- 


. Eesketh 


a curve  shape  of  the  fora  exp  P -(  V - L)  1 

L laax  J 

and  Schneider  find  /y\_  * 2 on  the  high  frequency  side,  but  /YV.  = 2. 5 
on  the  low  frequency  side.  The  latter  may  be  partly  a result  of 
anharmonic  forces. 

It  has  been  pointed  out  to  the  author  by  D.  L.  Dexter  that  the 
customary  use  of  Smakula’s  formula  to  determine  the  number  of  F centers 
in  a crystal,  or  their  oscillator  strength  is  incorrect  because  cf 
Smakula's  assumption  of  a Lorantzian  shape,  so  that  a considerable 
amount  of  experimental  data  must  be  reevaluated.  If  we  note  that  the 
absorption  constant  0^(  ) equals  (\Jq  ^ where 

A/c  is  the  number  of  centers,  the  intecrated  absorption  for  a line 
of  oscillator  strength  f can  from  (2. 5-2.7)  be  written: 

2TT  1a  e (3.30) 


/ha 


where  '0:  /V\.  = square  of  index  of  refraction.  Unfortunately 

the  usual  evaluation  of  experimental  data  estimates  the  integrated 
absorption,  not  by  integration,  but  merely  by  taking  a suitable 
produced  of  maximum  absorption  full  width  at  half 

power  W . For  a Lorentzian  shape 

I.5Tc(„,W  <3.31) 


whereas  for  a Gaussian  shape 


jot!  \.0Mm\Al  (3-3a) 
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Snakuia's  formula  is  equivalent  to  the  use  of  (3 *30)  and  (3*31)  combined 
with  the  assumption  (Eg/E)  2=.  ( /Tv.  + 2)  /9,  i.e.  the  assumption  that 
the  effective  field  E^  is  equal  to  the  Lorentz  field  E (^n/i)  P* 

On  the  basis  of  the  experimental  shape,  (3*30)  should  be  more  nearly 
combined  with  (3*32).  Furthermore,  the  F center  electron  is  sufficiently 
diffuse  that  the  effective  field  is  considerably  less  than  the  Lorentzian 
field.  A more  detailed  discussion  of  these  points  will  be  given  in  the 
Physical  Review  by  D.  L.  Dexter. 
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if.  Broadening:  Relation  of  Many  Body  and  One  Body  Viewpoints 


Although  the  one  parameter  treatment  discussed  above  is  an  over- 
simplification of  the  situation  for  a real  impurity  center,  it  is 
nevertheless  quite  feasible  to  use  a simple  configuration  diagram, 
for  the  purpose  of  deciding  whether  a semi-classical  treatment  is 
permissible  and  to  understand  the  observed  spectrum.  We  shall  try  to 
explain  in  terms  of  a many  body  viewpoint  why  for  most  problems  a one- 
I>arameter  analysis  is  likely  to  yield  qualitatively  correct  results. 

Our  essential  conclusion  will  be  that  if  the  modes  which  participate  in 
the  broadening  do  not  have  too  groat  a frequency  spread,  a one-parameter 
viewpoint  using  e suitable  mean  frequency  will  yield  fairly  accurate 
results. 

We  shall  start  by  assuming  that  the  configurational  energy 
difference  is  a function,  possibly  non-linear,  of  a single  parameter 
X where  the  latter  is  a linear  combination  of  the  displacements  or  the 


normal  coordinates  Q:  of  the  initial  state: 

O', 5- 

y' 


ih.l) 


where /V  is  the  total  number  of  normal  coordinates. 

If  the  semi -classical  approximation  is  valid,  we  proved  in  I 

that  the  one-parameter  spectrum  (3-6)  and  (3*7)  is  valid  (c.f.  1(5*7) 

2 

and  l(5*19))  providing  merely  that  x > Is  given  its  correct  rather 
than  its  one  parameter  value.  This  value  from  (4.1)  is 
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The  question  -/e  oay  ask  is  this:  is  there  a harDC;iiic  oscillator  with 

juenc 

K/ 


effective  mass  and  frequency  such  that  its  mean  oscillation 


(4.3) 


is  a good  approximation  to  (4.2)  at  all  temperatures?  Both  formulas  are 

linear  with  temperature  at  high  temperatures  and  have  a zero-point  value 

at  low  temperatures.  The  best  timt  we  can  do  is  to  choose  M and  cu 

6 e 

SC;  that  there  is  agreement  with  the  high  temperature  slope  and  the 
zero  point  value.  This  leads  to  the  choice 

, ivA/_\  ^ v-i\  /rA.^f. 

\'^  J/  io->  ij 


^ I « « 


“e  = r*  /Hj 


^4.4; 

(4.5) 

where  R = ^ /Jy  (}  .6) 

The  averaging  operation  ^v“  used  in  (4.4)  and  (4.5)  is  a weighted 


average  over  the  modes  using  Q.^  as  weight: 


/tv  ^^3=  2 


(‘>■7) 


It  is  clear  that  if  the  coefficients  that  measure  the  extent  to 


which  modes 


% 


participate  in  the  broadening  emphasize  a small  range 


of  frequencies,  then  m will  be  somewhere  within  the  range,  and  will 
be  essentially  M/R.  There  is  a certain  amount  of  genuine  arbitrariness 
in  Mg  through  the  choice  of  the  normalization  constant  R.  For  example 
if  X represents  the  s«un  of  six  radial  dlspi^rxements,  then  the  kinetic 
energy  would  by  6 M(x/.6)^  = M/6  and  R = 6.  If  x represents  the 
average  of  the  six  radial  displacements,  then  the  kinetic  energy  is 
6 M X and  R = 1/6.  In  any  case,  M/R  is  equivalent  to  the  intuitive 
choice  for  the  effective  mass.  The  presence  of  dispersion  among  the 
important  frequencies  can  only  increase  the  effective  mass  since 


/hr  fujV  ^ ^ 


Z 
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To  illuFtrate  the  above  remarks,  we  shall  consider  two  cases.  In 
the  first,  we  shall  let  x = , the  displacement  of  one  atom.  In 

this  case  all  the  are  equal.  For  the  acoustic  modes,  with  a 

Debye  spectrum  this  yields 


= (o)/2) 

(Mg/M)  = U/(3R) 


(4.8) 

(4.9) 


where  25  is  the  highest  acoustic  frequency.  If,  however,  we  were  to 
consider  a coordinate  x = xj^  - say  the  difference  between  the 

displacement  of  a K atom  on  each  side  of  the  T1  center  in  Williams  case, 
the  weighting  factor  associated  with  a phonon  of  propagation  constant 


would  be  1 1 - cos 


^ ■ ^-l)J 


which  takes  account  of  correlations 


in  the  motions  of  the  two  particles-  The  contributions  then  come 
largely  from  short  wave-length  acoustical  and  optical  modes.  (For  long 


wave  length  modes  of  either  type,  the  two  K atoms  move  iu  unlsou). 


\ 


only  acoustical  modes  were  considered,  we  would  find  ^ *66  (5,  and 

1.1  (M/R).  This  estimate  assumes  a Debye  spectrum  unperturbed  by 
the  center.  A stiffness  increase  near  the  center,  and  the  inclusion  of 
optical  modes  will  both  tend  to  raise  cu^.  Note  that  the  modes 
which  are  used  in  (3*29)  should  be  the  perturbed  modes  of  the  lattice 
in  the  presence  of  the  center.’  If  the  stiffness  change  near  the 
center  is  sufficiently  large,  a trapped  molecular  type  mode  will  occur, 
and  this  one  mode  will  make  a predominant  contribution  to  the  parameter 
X,  thus  insuring  the  validity  of  the  one-parameter  viewpoint. 

When  a fully  quantum  mechanical  calculation  is  necessary,  it  is 
more  difficult  to  assess  the  validity  of  a one-parameter  viewpoint, 
since  the  quantum-mechanical  many -body  problem  has  not  been  completely 


solved.  Three  cases  have  been  solved  however; 


(1)  The  energy  difference  ^ is  linear  in  the  normal  coordinates, 

29. 

and  the  latter  all  have  the  same  frequency  (Huang  and  Rhys  ) . 

(2)  The  energy  difference  ^ is  linear  in  the  normal  coordinates,  and 

18 

the  latter  have  a spread  in  frequencies  (Lax  ). 

(3)  The  energy  difference  contains  a diagonal  quadratic  form 

■i — Tl. 

in  addition  to  linear  terras 


(O'Rourke^®). 

The  pertinent  results  may  he  summarized  as  follows; 

(a)  The  Huang  Rhys  spectrum  for  many  modes  of  the  same  frequency  is 
identical  to  that  which  would  have  been  ohlalned  from  a eingip 
parameter  with  the  same  frequency  and  a suitably  chosen  linear 
energy  difference. 

(b)  The  Lax  spectrum  1(3.12)  reduces  to  the  Huang  and  Rhys  spectrum 
1(9 -l)  i'f  there  is  only  a small  dispersion  of  frequencies--so  that 
a one-parameter  description  is  possible. 

(c)  The  O'Rourke  spectrum  is  identical  in  shape  to  the  Lax  spectrum 

but  has  a shift  of  mean  energy  of  an  amount  . R^' 

This  shift  is  exactly  what  one  would  expect  from  the  moment 
formula  1(4.3)*  The  same  frequency  shift  can  be  obtained,  to  a 
good  approximation  in  a one-parameter  problem  by  having  a suitable 
stiffness  change.  In  a one  parameter  problem  however,  this 
frequency  shift  will  alwpys  be  accompanied  by  a change  of  shape. 

In  this  case,  it  is  not  the  one-parameter  result  which  is  wrong, 
but  rather  that  a quadratic  change  in  potential  energy  near  a 


I 


center  P)  X becomes  a quadratic  form 

containing  important  off-diagonal  terms  when  expressed  in  terras 
of  the  normal  coordinates. 


R.  C.  O'Rourke,  Phys.  Rev.  £1,  265  (1953)* 
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Pending  the  results  of  a nore  detailed  analysis  involving  off- 

diagonal  quadratic  terms,  it  is  probably  safe  to  say  that  the  criteria 

for  the  validity  of  a one  parameter  viewpoint  will  remain  valid: 

The  important  modes  that  contribute  must  hove  a small  dispersion  of 

frequencies.  If  this  is  the  case,  it  will  be  possible  to  choose  an 

effective  frequency  and  mass  to  fit  the  breadth  of  the  spectrum  over 

the  entire  temperature  range.  If  appreciable  skewness  is  present, 

however,  it  may  be  difficult  to  fit  over  an  entire  temperature  range 

by  choosing  one  parameter  B,  unless  some  restrictinns  are  pieced  or. 

the  second  order  energy  change  f\i~'  ^ , Bej. 

As  an  interesting  example  of  the  relationship  between  the  one 

parameter  and  the  many  body  viewpoint  ve  shall  discuss  the  case  of 

20 

Impurity  levels  in  semi-conductors.  We  note  that  for  S or on  doped 
Silicon  there  are  about  500  atcjms  within  a sphere  whose  radius  CX 
is  the  first  Bohr  radius  of  the  trapped  electron  in  its  Is  state. 

The  broadening  is  produced  by  the  vibration  of  these  atoms  through 
the  electron  lattice  interaction.  We  may  describe  the  latter  by  the 
deformation  potential.  ^ 

E,  di\)SRcr  1 E,  N (u.io) 


so  that  the  energy  shift  associated  with  an  electron  in  state 
is  given  approximately  by 


AE -( £, N (it.u) 


where 


= [j  f (xa-C (4.12) 
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and  the  explicit  form  for  A/^t“ y is  based  on  a hytirogenic  wave  function 
of  Bohr  radius  CX.  • The  normal  coordinate  associated  with 

phonons  of  propagation  constant  . 

To  connect  with  a one  p»araineter  viewpoint,  we  simply  introduce 
for  our  parameter  x,  the  energy  shift  itself.  The  constant  O-  ^ , 
now  written  O-t;;  obey 

|ar\’^=  ^ (It. 13) 


:ffective  frequency  assuming  a Debye  spectrum  can  now  be  obtained 


from  (t.h)  and  (4.11); 




1 (o  \T 
3tt  <x 


where  O'"  is  the  Velocity  of  sound.  The  Integrals  in  (4.l4) 

were  evaluated  by  extending  the  limits  to  infinity  since  -L  Ou^ 

1.  rr(\o^ 

7.7. 

The  (squared)  broadening  at  absolute  zero  may  be  obtained  by 
setting  ^ = O (^.2); 

- z (‘••15) 


Mvcl 


The  ratio  of  the  squared  broadening  ^ to  its  absolute  zero 


value  can  froTn  (4.2)  be  written: 


C d.  -* 

pfr)-  \ + 24-  \ x^(i  +x^r  ^ V Tj  ■ ■ Q cix 
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where  x = ^ CX.  ^ and  ^ = 2 KT  / 'fb  .An  approxiiaate 

formula  for  R(T)  using  the  single  effective  frequency  (o^  is  simply 

R,(t)  = 1 + 2 =■  Cote,  (k  ^ 'oJe  /^t) 

We  note  that  “^CD^/k  = (8/3«)A.  A comparison  of  the  "exact"  and 
approximate  R(T)  from  (4.l6)  and  (4.17)  is  shown  in  Figure  8.  They 
both  have  the  same  zero  point  value  and  high  temperature  slope.  Near 

T = o R(tW  1 H-  O.'l  (z  — T/a  r 

vhe«a.  ^ 1 +-Z  (-  + A/3fYT') 

The  exponential  decrease  in  R (T)  is  characteristic  of  the  Einstein 

approximation  we  have  made,  whereas  the  power  law  decrease  in  R(T)  is 
characteristic  of  the  contributions  of  the  low  frequency  phonons— as 
in  the  Debye  theory  of  specific  heats. 

To  compare  these  results  with  experiment^  we  note  that  the  zero 
point  line  width  obtained  from  (4.15)  is  of  the  order  of  0.004  e.v. 
whereas  the  experimental  width  at  liquid  helium  temperatures  is  about 
0.001  e.v.  (see  figure  1).  The  theoretical  line  width  may  be  an 
overestimate  because  of  the  use  of  the  Born-Oppenheimer  approximation. 
It  may  also  be  an  overestimate  because  the  deformation  constant  has 
been  overestimated.  The  reason  for  the  latter  possibility  is  that  E^ 
is  determined  by  comparing  an  experimental  mobility  with  a theoretical 
mobility  based  on  thermal  scattering— neglecting  interband  transitions 
which  are  undoubtedly  significant  in  p type  silicon. 

To  see  whether  we  have  assumed  the  correct  basic  mechanism  for 
broadening,  namely  interaction  of  the  electron  with  acoustic  modes, 
we  can  examine  the  way  in  which  broadening  changes  with  temperature. 
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An  examination  of  (4.l6)  indicates  that  the  modes  of  importance  have 
an  energy  A where  A is  about  90“K  in  silicon.  These  are  modes  of 
quite  low  energy— they  have  wave-lengths  whose  size  is  of  the  order  of 
the  Bohr  radius  of  the  trapped  hole.  Does  experiment  bear  this  out? 
Figure  8 shows  that  R(T)  is  of  the  order  of  2 at  T = A,  i,e.  an 
appreciable  broadening  should  occur  theoretically  for  temperatures  as 
low  as  90®K.  An  examination  of  Figure  1 shows  that  a broadening  of 
just  the  right  order  of  magnitude^^^  occurs  between  nitrogen  and 
helium  temperatures.  This  substantiates  the  importance  of  low  energy 
modes  in  the  broadening.  Since  our  calculated  nhsolute  broadening  is 
of  the  right  order  of  magnitude,  our  basic  picture  of  the  phenomenon 
must  be  correct. 


5.  Thennal  Ionization  and  Capture 


As  is  undoubtedly  discussed  in  other  papers  of  this  conference,  the 
sensitivity  of  a photoconductor  increases  with  the  lifetime  of  the 
electrons  in  the  conduction  band,  but  its  rapidity  of  response  decreases. 
(Corresponding  remarks  can  be  made  about  holes  in  the  valence  band  but 
will  henceforth  be  omitted) . The  photoconductive  lifetime  is  determined 
by  a competition  between  the  rate  of  therral  ionization  and  the  rate  of 
recombination.  For  example,  if  we  are  dealing  with  a crystal  containing 
two  types  of  impurity  centers:  |\^  ^ donors/cm^  and  N|q_  acceptors/cm^. 


witn 


K I 


(\/q  > , then  we  will  liuVe  ii  ptictoccri'^iict ivit jT  witl? 


a lifetime 


(5.1) 


where  B is  the  recombination  rate  and  K = A/'a  is  the  ratio  of  ioaizutiou 
to  recombination  rates,  or  the  equilibrium  constant: 

K,  - (2.TT  E Mt)  '5-2) 


where  Eg  is  the  (ground  state)  ionization  energy  of  the  center,  assumed 
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to  have  only  one  state,  and  m is  the  effective  mass.  If  the  electron 
energy  has  several  minima  within  the  Brillouiu  zone,  and  an  ellipsoidal 


The  K of  (5.2)  must  be  divided  by  the  factor 

{1  + ^L( - EJ  M r]  I 

where  the  states  excited  states  and 

corresponding  degeneracy  factor.  See  Burstein,  Bell,  Davisson 
and  Lax,  J.  Phys.  Chem.  2Jf  (1953)  And  K.  S.  Shifrin, 
Teknicheskoi  Zhurniel  Fiski  l4,  ^3  (19^4). 
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energy  surface  with  masses  m2  and  m^,  then  ve  must  make  the 
replacement 

* l/p 

(m  ) > (No.  of  valleys)  (mj^m^m^)  (5.3) 

since  the  latter  expression  takes  proper  account  of  the  density  of 

states  near  the  mlnlar’m  band  energy. 

At  high  temperatures,  we  will  have  K >>  N^  and 


, (X  JioK^) 

the  Ionization  rate,  whereas  at  low  temperatures  we  will  have 


^ 6K=  A 


(5.i<> 


K <<  N.  and 

H 


4 ^ I 

C::  t!)  in/  A 


^ Y"  _ i\  r>\ 

\ \ iuyrr\OJUL ) 


) . . , 15.5) 

the  recombination  rate— since  the  number  of  Ionized  donors  at  low 


temperatures  will  be  approximately  equal  to  N the  number  of  acceptors 
available  to  remove  electrons  from  the  donor  levels.  In  germanium  with 

^ t. 

— lO'*’*^  acceptors/cm^,  the  temperature  dividing  the  regions  of 

validity  between  (5«^)  and  (5*5)  la  about  25 "K. 

As  we  shall  see  In  the  subsequent  discussion,  one  of  the  dlfflcultlee 

that  faces  a basic  calculation  of  Ionization  rates  Is  that  fairly  large 

values  must  be  obtained  In  order  to  agree  with  experiment.  For  example, 

14  3 

germanium  with  about  10  acceptoi-s/cm  has  a photoconductlve  lifetime 

_ -5  ^0 

^ of  less  than  10  seconds  even  at  liquid  helium  temperatures. 

This  linplies  that  B ^ 10  ^ cm^/sec,  or  since  K lO”^,  A > 10  ^^/sec. 

Perhaps  a more  understandable  way  to  put  this  statement  Is  to  note  that 

at  4*K,thermal  velocities  are  of  the  order  3 * 10^  cm/sec  so  that  the 

-16  2 

center  has  an  effective  cross-section  of  the  order  of  3 10  cm  or 

larger  for  thermal  capture. 
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Bursteln,  Davisson,  Bell,  Turner  and  Llpson,  Phys.  Rev.  ^5  (195^)- 
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As  a general  experinental  argument  for  high  radiationless 

transition  probabilities,  we  note  that  most  electron  traps  do  not 

luminesce.  In  particular,  a careful  search  fcr  F center  luminescence 

4l 

has  been  made  by  Klick  whose  essential  conclusion  is  that  radiative 
transitions  are  never  more  than  3 per  cent  of  the  non-radiative 
transitions. 

One  of  the  first  attempts  to  treat  thermal  ionization  from 

electrtns  trapped  in  semi-conductors  was  made  by  Goodman,  Lawson  and 
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Bchiff . They  acsumed  that,  the  essential  interaction  as  far  as  ioniza- 

tion is  concerned  is  the  interaction  between  the  trs-pped  electron  and  the 
vibrating  core  atom; 

\/  =1 ^ <5.6) 

Vc.Ls,  Kla-R\  KA  K 


where  R is  the  nuclear  displacement.  Their  unperturbed  wave -functions 
were  simple  ("Hartree")  products  cf  an  electronic -wave  function  times 
a vibrational  wave  function,  and  the  transition  is  caused  by  \J^  ^ ^ 
as  a perturbation.  We  shall  refer  to  this  as  the  Hartree  approximation. 

Using  only  terms  linear  in  R,  only  one-phonon  processes  are 
allowed  in  the  Hartree  approximation.  Since  the  ionization  energies  of 
hydrogen-like  traps  in  Ge  or  Si  are  less  than  the  corresponding  Debye 
energies,  acoustic  phonons  are  available  to  make  a one-phonon  transition 
possible.  The  Goodman-Lawson-Schiff  answer,  for  the  ionization 
probability  turns  out  too  small  by  about  two  orders  of  magnitude. 


C.  C.  Klick,  Phys.  Rev.  Q4,  1541  (1954). 

Goodman,  Lawson  and  Schiff,  Phys.  Rev.  191  (19^7)  • 
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Further  iiscuBsion  of  the  G.L.S.  interaction  will  be  given  later. 

Goodman,  Lawson,  and  Schiff  also  made  an  estimate  of  the 
probability  for  multi-phonon  processes  using  terms  of  order  in  the 

^+3  k2 

potential  and  found  their  contribution  to  be  small.  Seitz  in 
unpublished  work  treated  the  same  interaction  \l ^ in  a Born- 

Oppenheiner  approximation,  arriving  at  results,  according  to  Goodman, 
Lawson  and  Schiff  of  the  sane  order  as  their  own.  An  explanation  for 
this  agreement  will  be  given  later. 


During  LVic  ouiamcr  Ox  195^> 


Xil  i^XC<XUXU^  OiAC  kj  X Li  X 


4 *y**>»i'**4  ^ 


levels  in  silicon  with  the  deformation  potential  interaction,  I also 
estimated  the  thermal  ionization  probability  with  the  same  interaction 
using  the  Hartree  approximation.  Although  the  broadening  cane  out  one 
order  of  magnitude  too  high,  with  the  experimental  constants  then 
available,  the  thermal  ionization  was  an  order  of  magnitude  lew. 

Because  of  this  puzzling  discrepancy,  publication  of  this  work  was 
postponed,  and  both  problems  were  reexamined.  The  thermal  ionization 
problem  was  treated  in  the  Born-Oppenheimer  approximation  and  the  many 


The  problem  of  non-radiative  transitions  between  bands  had  been 
discussed  .ourlier  by  F.  Moglich  and  R.  Rempe,  Physikalische  Zeitschrif 
4l,  236  {l9kQ);  Zelts.  fur  Phys.  li;?,  707  (19^0).  They  used  a 
Hartree  viewpoint  and  suggested  that  the  many-phonon  processes 
required  here  by  energy  conservation  wore  obtained  from  terms  In 
the  electron  lattice  interaction  of  high  order  in  the  nuclear 
displacements.  Only  qualitative  estimates  were  made  for  transition 
rates. 
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body  suras  were  performed  following  raethods  used  in  the  broadening 
l8 

analysis.  But  the  answer  did  not  change  appreciably.  The  disagreeraent 
for  the  brcjadening  calculation  was  reduced  in  the  raeantirae  partly  by 
taking  into  account  "motional  narr^vl ng"  and  partly  by  changed  . 

44 

experimental  constants  due  to  Improved  nobility  measurements. 

L5 

Leurgans,  in  his  thesis,  also  considered  the  thermal  ionization 

problem  with  the  deformation  potential  and  the  Hartree  approximation. 

46 

His  results  confirmed  the  discrepancy  we  found.  Yafet  reconsidered 
xhe  problem  in  the  Born-Oppenheimer  approximation  and  found  that  multi- 
phonon contributions  were  negligible  sc  that  the  answer  Leurgans  and 
I had  obtained  using  the  Hartree  approximation  would  nc.t  be  appreciably 
increased  by  raulti-phonoin  contributions.  An  appreciable  increase  was 
found,  however,  when  Coulomb  wave -functions  rather  than  plane  wavea  were 
used  for  the  final  state.  A more  detailed  discussion  will  be  given 

later  in  this  section. 

47 

Kubo  suggested  that  perhaps  quadratic  terras  in  the  interaction 
energy  are  iraportant--so  that  the  vibration  frequencies  in  the  excited 
and  ground  electronic  states  were  different.  He  used  the  Goodman, 

Lawson  and  Schiff  interaction,  the  Born-Oppenheimer  approxiraation  and 
treated  the  central  atom  as  a simple  harmonic  oscillator.  He  more  or 
less  arbitrarily  chose  the  change  in  frequency  of  the  harmonic  oscillator 
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M.  B.  Prince,  Phys.  Rev.  93,  1204  (1954). 

P.  J.  Leurgans,  Thesis  1952  University  of  Illinois. 
Y.  Yafet,  Private  communication. 

R.  Kubo,  Phys.  Rev,  80,  929  (1952). 
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because  of  a change  in  electronic  state  to  be  10  per  cent.  Because  of 
the  arbitrariness  of  Kubo's  model,  it  is  difficult  to  assess  the 
applicability  of  his  results,  however,  there  is  no  question  but  that 
any  significant  change  in  vibrational  frequencies  will  appreciably 
increase  the  rate  of  themal  ionization. 

In  addition  to  the  case  of  electrons  trapped  at  impurity  levels 
in  semi-conductors,  considerable  effort  has  been  spent  on  radiationless 
transitions  in  F centers.  From  the  formal  point  of  view,  the  first 
major  contribution  was  made  by  Huang  and  Rhys  who  made  a many  body 
calculation  in  the  Born-Oppenheimer  approximation  assuming  an  interaction 
energy  linear  in  the  normal  coordinates.  The  major  limitation  of  their 
calculational  technique— as  in  their  treatment  of  broadening--is  that 
all  the  modes  must  have  one  frequency.  (This  limitation  can  paally  be 
by-passed  by  using  techniques  developed  in  the  author's  analysis  of 
the  Franck-Condon  principle). 

The  rate  for  radiationless  transitions  from  the  2p  to  the  Is  state 
calculated  by  Huang  and  Rhys^  was  so  low  in  comparison  to  radiative 
transitions,  (lO"^^  in  KBr  at  30“K)  that  luminescence  should  be  easily 
observed  at  most  temperatures  for  all  the  alkali  halides.  This  is 
definitely  not  the  case,  however. 

he 

Meyer  suggested  that  a possible  explanation  for  the  low 
radiationless  transition  rate  was  the  neglect  of  quadratic  terms  in  the 
energy  difference  between  the  two  states.  He  introduced  therefore,  a 
diagonal  quadratic  form  in  addition  to  the  Huang-Rhys  interaction  energy 
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H.  J.  G.  Meyer,  Physica  l8l  L (1954). 


and  recalculated  the  transition  rate.  The  diagonal  coefficients  are 
somewhat  arbitrary,  how’ever>  and  he  chose  then  so  as  to  agree  with  the 
ohsei-vcd  temperature  shift  in  the  optical  absorption  peak.  He  then 
found  a sufficiently  high  radiationless  rate  to  explain  why  luminescence 
h>is  been  difficult  to  observe.  (Badiation  less /radiative  Cr  30  KBr 
at  20°K). 

In  a subsequent  unpublished  manuscript,  however,  Meyer  comes  to 
entirely  different  conclusions: 

(1)  He  finds  it  permissible  to  neglect  quadratic  terms  in  the  energy 
difference— having  first  estimated  their  size  by  comparison  with 
experiment. 

(2)  He  finds  that  the-  ratio  of  radiationless  tc  radiative  transition 

probabilities  from  the  excited  to  the  ground  stats  of  the  F center 
-7 

is  now  10  in  KBr  at  0®K.— so  that  luminescence  should  be  observable. 
The  difference  between  this  result  and  his  previous  one  is 
apparently  due  to  different  ways  of  determining  empirically  the  Huang- 
Rhys  parameter  S which  measures  the  strength  of  the  linear  terms  in  the 
interaction  energy.  If  H is  the  full  width  at  half -power,  his  new 
method  of  determining  S is  to  use  the  relationship: 

which  is  found  to  hold  empirically.  The  formula  previously  used,  which 

had  a thccretical  justification  was  identical  to  (5*7)  with  C = 0. 

It  was  not  clear  to  me  why  a moderate  change  in  S would  change  his 

8 

answer  for  radiationless  transitions  by  a factor  10  . However,  a 
detailed  examination  of  his  formula  in  the  limit  T = 0°K  revealed  a 
factor  where  is  the  number  of  phonons  involved  in  the 
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transltion--which  is  of  the  order  70  in  KBr--sc  that  a 20  per  cent 

g 

error  in  S can  yield  a factor  10  in  the  transition  probability;  In 

view  of  the  sensitivity  of  Meyer's  result,  conclusions  about  the 

feasibility  of  observing  luminescence  must  be  regarded  as  provisional. 

An  attempt  to  generalize  the  Huang-Rhys  method  to  the  case  when  a 

k9 

distribution  of  frequencies  are  present  has  been  made  by  Tewordt  for 

application  to  semi-conductors.  Tewordt 's  approach,  however,  parallels 

Huang  and  Rhys  closely,  and  he  is  consequently  unable  to  carry  out 

explicitly  the  sums  over  the  vibrational  states.  He  was  apparently 

unaware  that  this  portion  of  the  problem  had  already  been  solved. 

50 

Pekar  has  written  an  interesting  paper  on  the  diffusion  of 
electrons  toward  a trap  and  their  eventual  capture.  His  analysis  uses 
classical  statistics  and  the  diffusion  equation  outside  a somewhat 
arbitrary  radius  at  which  capture  is  supposed  to  take  place  at 

a rate  ^ where  /w(/lo^  density  at  . His 

3 

final  result  is  an  effective  capture  rate  cm  /sec  given  by: 


(5.8) 


1 + 


^ p ’/Ao 

^ (4^D)'V^  [-e  Vf«.)/RT] 


^ L.  Tewordt,  Zeits  fux'  Physik  1^,  6ok  (195^). 

S.  I.  Pekar,  Abh.  Sow.  Phys.  1,  47  (1951)  translated  from 
Z.  f.  exp.  u.  theor.  Phys.  267  (1950). 
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If  the  second  terra  in  the  denominator  is  large  compared  to  the  first 
(capture  rate  large  compared  to  diffusion  rate),  the  process  vill  be 
diffusion  limited  with  a rate  given  by: 

d ( k)]'"  (5.9) 


The  diffusion  constant  D can  be  calculated  from  the  mobility  Jj.  using 

51 


the  Einstein  relation 


' c-  .1  -.s  \ 


The  potential  can  be  taken  to  be  eV(  /\.  ) = "(e  / Kv  i\/  ) where  is 
the  dielectric  constant.  At  sufficiently  low  temperatures  that 

integral  in  (5*9)  becomes 

insensitive  to  its  upper  limit  which  may  be  set  equal  to  infinity.  Thus 
if  -pjT  « Eg  = the  ionization  energy  of  the  center,  we  find  that 

Bd  - H-lTy0.e/t^  (5-11) 

The  diffusion  rates  thus  obtained  are  quite  high.  At  liquid  heliura 

5 

temperatures,  using  an  experimental  mobility  of  the  order  of  5 x 10  , 

*•2. 

we  find  0,  J(  Cm  6 x 10  cm'^/sec , whereas  the  calculated  quantum 
capture  rate  is  or  the  order  of  0 ^ 10  cm  /sec  in  germanium  at  the 
same  temperature.  Thus  the  diffusion  is  very  rapid,  and  the  limiting 
factor  is  the  capture  p?.*ocess  itself. 


For  experimental  verification  see  Transistor  Teachers  Summer  School, 
Phys.  Rev.  I368  (1952). 

52  E.  Conwell,  Proc.  I.  R.  E.  40,  132?  (1952). 
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It  is  of  some  interest  to  note  that  if  we  define 

g)  = jS  e \/('u^  I (5.12) 

to  he  the  capture  rate  adjusted  for  the  increase  in  density  near  the 
oricin,  Pekars  equation  can  be  rewritten  in  the  fora 

J_  _ _L_  J L '5.13) 


1)0 


B 6a 


indicatinc  that  the  total  time  for  trapping  is  the  sum  of  the  diffusion 
tiae  and  the  adjusted  capture  time.  Since  the  diffusion  time  is,  by 
uui-  estiaaten,  'sually  negligible,  we  shall  be  concerned  primarily  with 
the  adjusted  capture . time . 

The  adjustment  (5*12)  is  somewhat  sensitive  to  the  choice  of  j 

but  the  latter  can  be  taken  to  be  of  the  order  of  a Bohr  radius  so  that 
Cr  |6)  exp(Eg/kT).  Since  we  have  assumed  kT  is  siuall  compared 
to  the  ionization  energy  Eg,  an  enormous  increase  in  capture  would 
occur,  if  this  result  were  taken  literally.  However,  the  electron 
density  cannot  suffer  such  a marked  increase  near  the  center  because  of 
the  uncertainty  principle. 

It  is  therefore  necessary  to  take  account,  quantum  mechanically, 
of  the  increase  in  density  near  the  center.  G'Jimmel  and  I have  done  this 
by  using  Coulomb  continuum  wave  functions  instead  of  plane  waves  to 

53 

describe  the  incoming  electrons.  Roughl.y  speaking  the  increase  in 
density  near  the  center  is  the  Sommerfold  factor.  A detailed 
evaluation  of  the  matrix  element  for  the  transition  yields  an  additional 
factor  2,  so  that  the  appi-opriate  correction  in  the  squared  matrix 
element  for  the  use  of  coulomb  rather  than  plane  wav-js  is 


H.  Gummel  and  M.  Lax,  Phys.  Rev.,  to  be  submitted. 
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(8«/  ^ ■ exp(-2n/^  )"^  ^ where  ka,  k is  the  propagation 

constant  of  the  incident  electron,  and  a is  the  Bohr  radius  of  the  trapped 
state . 

A rough  Indication  of  the  Coulomb  correction  at  low  temperatures 
can  be  obtained  by  noting  that  the  important  k's  in  the  equilibrium 

distribution  are  of  the  order  ( k /2m*^  ) ^2=1-  k^  T or 

l/2  l/2 

ka  r\^  (kpT/Eg)  ' . Thus  the  correction  factor  is  roughly  Gjr(Eg/koT) 

which  can  be  of  the  order  of  200  at  liquid  helium  temperatures  (for 

Eg  .0125  e.v.  in  germanium).  At  high  temperatures,  the  Boltzmann 

factor  puts  only  a weak  limit  on  k,  but  the  matrix  element  for  the 

deformation  potential  interaction  contains  a factor  |^1  - (ka)^ J 

which  means  the  important  values  of  ka  are  smaller  than  say  1/2.  Even 

at  high  temperatures,  then,  the  correction  factor  can  be  as  large  as  50. 

To  illustrate  these  results  in  more  detail,  we  note  that  the  total 
capture  rate  0)  is  a Boltzmann  average  over  the  capture  rates 
from  a given  state  k; 


on 


(5.14) 


For  simplicity,  we  assume  spherical  energy  surfaces,  so  that  the 
normalized  distribution  function  is  given  by: 

The  transition  probabilities  can  be  written  in  the  form: 


"t: 


;i.TrK  Into  O' 


{ I V rr*)  O 


(5.15) 


(5.26) 
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where  A is  the  voiune  of  a unit  cell,  is  the  raass  of  one  atom, 

f is  the  propagation  constant  of  a phonon,  ^ o)  = cn(  "t'  ) is 
its  energy,  and  VJ"  = dco/d'^is  the  corresponding  velocity  of  sound.  We 
must  insert  for  and  m a value  consistent  with  conser/ation  of  energy: 

Ato  (5.17) 

If  we  let  cu(  ) = "f\  “o  “ then  for  a hydrogenic  impurity 

(5*17)  is  equivalent  to  the  conditions 

\ U3  - U)J  1 ) (5.18) 

For  the  case  in  which  the  perturbing  potential  is  the  deformation 
potential  VC  have  * . 

Q -=  i £,  I ( e' " (5.19) 

The  matrix  element  of  exp(  between  the  free  state  k and  the 

bound  state  a has  been  calculated  exactly.  Because  AT  is  large  compared 
to  a“'^,  the  reciprocal  Bohr  radius,  however,  the  resulting  ejepression  can 
be  greatly  simplified.  If  the  free  state  is  described  by  a plane  wave 


we  get 


"T.vJ. 


a ' “ \ 


(5.20) 


If  the  free  state  Is  described  by  a Coulomb  wave  we  get 


(5.21) 


Thus  justifying  the  previous  remarks  about  the  importance  of  the 
Coulomb  correction. 

The  corresponding  formulas  for  D using  the  Goodman,  Lawson,  Schiff 
interaction  are: 
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(5.22) 


n ^ 3.2TLa 


V 


n - ^TT  g ii/6pr~  (%-)  W~'v^l 

4»j , ' V u.  (i  t if  )ii-  s>4^  f- 


(5.25) 


The  four  values  of  D can  be  inserted  into  (5*16)  and  (5*1^)  to  obtain 
expressions  for  the  total  capture  rate.  For  simplicity,  ve  list  the 
results  for  the  temperature  range  k^jT  Eg — the  range  of  greatest 

interest. 


lUE.n  1 

MEjiro- 


(5.24) 


(5.25) 


0 
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P.W, 


g(^T)Q(t:,a) 

3 vro- 


(5.26) 


fc4-'3TTetM^(r.4' 

M'vo^ 


(5.27) 
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At  liquid  helium  temperatures  In  germanium,  the  capture  rate  using 

..<5  -3 

the  deformation  potential  and  Coulomb  waves  is  about  10  ' cm'^/sec. 

The  use  of  plane  waves  yields  a result  smaller  by  a factor  of  about  175* 
The  Goodman-Lawson-Schitf  interaction  with  plane  waves  gives  a result 
smaller  by  a factor  5000;  with  Coulomb  waves  the  result  is  smaller  by 
a factor  23.  Thus  the  primary  mechanisui  at  low  temperatures  seems  to 
be  the  usual  electron  lattice  interaction--and  Coulomb  corrections  are 
Important . 

A B of  10  ^ cm^/sec  with  10^^  compensating  impurities  would  lead 


to  a ‘*^hotoccnduct  ive 


life-time 


of 


- 
I O 


(5.28) 


On  the  other  hand,  Bursteia  et  al,' 

-5 


Uo 


find  that  for  Za  doped  germanium 


10  ■ sec  is  an  upper  limit  for  the  life-time— so  that  the  theoretical 
capture  rate  may  be  all  right,  but  is  possibly  too  small. 

For  silicon,  the  theoretical  capture  rate  is  10  ca^/eec.  The 
large  increase  compared  to  germanium  is  produced  by  a smaller  value  of 
t;a.  because  of  the  higher  velocity  of  sound  in  silicon.  Assuming 
the  number  of  compensating  impurities  in  fairly  pure  silicon  to  be  t>f 
the  order  of  10^^/cm^,  the  photoconduct ive  life-time  is  given  by 

'tr  ^ (10  10  )"  iirr  lo"  sec  which  is  well  within  the  experimental 

-4  5^+ 

limit  ^ < 10  sec  set  by  Burstein. 


5^ 


E.  Burstein,  J.  J.  Oborly  and  J.  W.  Davisson,  Phys.  Rev.  331  (1953)' 
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Infrared  absorption  of  a single  crystal  of  silicon  doped  with  boron  (E.  Bur  stein, 
E.  E.  Bell,  and  B.  W.  Hen  vis,  unpublished).  Further  discussion  of  this  data  may 
be  found  in  J.  Phys,  Chem.  57,  849  (1953),  (Burstein,  Bell,  IDavisson,  and  Lax). 


Infrared  absorption  of  a single  crystal  of  silicon  doped  with  boron  (E.  Burstein, 
E.  E.  Bell,  and  B.  W.  Heiivis,  unpublished).  Further  discussion  of  this  data  may 
be  found  in  J.  Phys.  Chein.  57,  849  (1953),  (Burstein,  Bell,  Davisson,  and  Lax). 
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Figure  3-  Configuration  coordinate  diagram  of  the  activator  system  for 
KC1:T1  including  the  effect  of  angular  dependence  of  the  Tl'*’ 
charge  density  and  the  perturbation  of  the  states  of  T1  by 
the  crystal  field.  (F.  E.  Williams,  J.  Phys.  Chem.  ?80 

(1953)). 
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Figure  4.  Theoretical  and  experimental  absorption  and  emission  spectra 
of  KC1:T1.  (F.  E.  Williams,  J.  Phye . Chem.  T80  (1953))* 
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Figure  5-  The  influence  of  temperature  on  the  absorption  spectrum  of 
color  centers  in  KBr.  (E.  Mollwo,  Zeits.  fur  Physik  8p,  56 


The  theoretical  ratio  R(T)  of  the  squared  broadening  at  the  temperature  T 
to  its  zero  point  value,  is  plotted  as  a function  of  the  dimensionless  variable 
x/d  . The  temperature  A = 2fjv/kpaw  94°K  in  silicon,  where  v is  velocity  of 
sound,  a is  the  Bohr  radius  in  the  crystal,  and  kg  is  the  Boltzman  constant. 
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